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Abstract. A factorization of a permutation into transpositions is called “primitive” if its factors are weakly ordered. 
We discuss the problem of enumerating primitive factorizations of peimutations, and its place in the hierarchy of pre¬ 
viously studied factorization problems. Several formulas enumerating minimal primitive and possibly non-minimal 
primitive factorizations are presented, and interesting connections with Jucys-Murphy elements, symmetric group 
characters, and matrix models are described. 

Resume. Une factorisation en transpositions d'une permutation est dite “primitive” si ses facteurs sont ordonnes. 
Nous discutons du probleme de remuneration des factorisations primitives de permutations, et de sa place dans 
la hierarchie des problemes de factorisation precedemment etudies. Nous presentons plusieurs foimules enumerant 
certaines classes de factorisations primitives, et nous soulignons des connexions interessantes avec les elements Jucys- 
Murphy, les caracteres des groupes symetriques, et les modeles de matrices. 

Keywords: Primitive factorizations, Jucys-Murphy elements, matrix integrals. 


1 Introduction 

1.1 Polynomial integrals on unitary groups 

Let U{N) denote the group of N x N complex unitary matrices U = [uij\i<i,j<N- By a polynomial 
function on U(N) we mean a function of the form 

p(u)= E (i) 

m,ri >0 /,«/,/' ,J' 

where 

I (? 1 , ■ • • , i m ) I (* 1 > • • • ! * 71 ) 

J = (jl> • • • ; jm) J = Ol) • • ■ ! Jn) 

are multi-indices, 

UijUpji Ui 1 j 1 . . . 
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( 2 ) 

( 3 ) 
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is the corresponding monomial in matrix entries, and only finitely many of the coefficients c(/, J, I', J') £ 
C are non-zero. A polynomial integral over U(N ) is the integral of a polynomial function on U(N) 
against the normalized Haar measure. 

The computation of polynomial integrals over U(N) is of interest from many points of view, including 
mathematical physics (nuclear physics, lattice gauge theory, quantum transport and quantum informa¬ 
tion), random matrix theory (matrix models, asymptotic freeness of random matrices), number theory 
(stochastic models of the Riemann zeta function), and algebraic combinatorics (integral representations 
of structure constants in the ring of symmetric functions), see m for references to the large body of 
literature on matrix integrals of this type. Nevertheless, the evaluation of such integrals is a problem of 
substantial complexity that is not yet fully understood. 

We wish to develop a general theory of polynomial integrals over U(N). By linearity of the integral, 
we have 

f p(U)dU= V V f U I jU I .j.dU , (4) 

JU(N) m,n>0I,J,I',J' JU{N ) 

so we consider the problem of evaluating monomial integrals 


/ UuUpj.dU. (5) 

JU(N) 

Monomial integrals are already of great interest in mathematical physics, see e.g. 0 . An easy argument 
involving the invariance of Haar measure shows that 0 can be non-zero only for m = n (i.e. the multi¬ 
indices I , J are of the same length as the multi-indices /', J'). Furthermore, when m = n < N (i.e. the 
degree of the monomial to be integrated is at most the dimension of the matrices being integrated over), 
the integral ([5} can be decomposed into a double sum over the symmetric group S(n) of the form 



UijUpj'dU 


£ [I = a{I')}[J = r{J')} W, 

(cr,r)G-S'(n) xS(n) 


( 6 ) 


This integration formula has two ingredients: a combinatorial “Wick-like” rule — sum over pairs of per¬ 
mutations (cr, t) such that a maps the multi-index /' to the multi-index / and r maps the multi-index J' to 
the multi-index J — together with a certain “weight” W CTr associated to each admissible pair of permuta¬ 
tions. These weights have a remarkable combinatorial interpretation as generating functions enumerating 
certain factorizations in the symmetric group; the resulting connections with algebraic combinatorics are 
the focus of this extended abstract prepared by the authors for FPSAC 2010. 

1.2 Primitive factorizations and Weingarten numbers 

Let S( oo) denote the group of finitary permutations of the natural numbers { 1 , 2, 3,... }, with S(n) < 
S(o o) the subgroup of permutations of [1, n] = {1,..., n}. An ordered sequence of transpositions 


(^1 fl )(^2 ^ 2 ) • ■ • ffc )) < L , 
is said to be a.factorization of n £ 8 ( 00 ) if 


(7) 


7T = (Sl ft) O (s 2 t 2 ) O • • • O (S k t k ). 


( 8 ) 
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ti < t 2 < ■ ■ ■ < t k 


(9) 


hold in (| 8 j. Consider the quantities 


hk.n{ti) = ^{factorizations of n into k transpositions from S(n)} 

Wk,ir(n) = ^{primitive factorizations of n into k transpositions from S(n)}. 


( 10 ) 


The numbers hk^(n) are known as (disconnected) Hurwitz numbers, and are of much interest in enumer- 
ative geometry, see e.g. ma. We will call the numbers u>fc )7r (n) Weingarten numbers, see El do) for the 
origin of this name. The primitive factorizations counted by Weingarten numbers have previously been 
considered by combinatorialists, both in relation to the enumeration of chains in noncrossing partition lat¬ 
tices III 01 and for their own sake 0. Our approach to polynomial integrals over unitary groups is based 
on the remarkable fact that the weights appearing in the integration formula { 6 } are generating functions 
for Weingarten numbers. 

Theorem 1 ( iflOllTIll ) For any n < N and n £ S(ir) we have 



where n = cr o t 


2 Jucys-Murphy elements 

2.1 Centrality 


Let C ;/ C S( oo) denote the conjugacy class of permutations of reduced cycle type p (p is a Young 
diagram). For instance, C( 0 ) is the class of the identity permutation, C(i) is the class of transpositions, 
and more generally Ct r ) is the class of (r + l)-cycles. Note that \p\ is the minimal length of a factorization 
of 7r into transpositions. The conjugacy classes of S(n) are C^(n) := C D S(n). Let Z{n) denote the 
centre of the group algebra C[5(n)]. Then {C^(ri')} is the canonical basis of Z(n ), where C^in) is 
identified with the formal sum of its elements, so Z(n) is referred to as the class algebra of S(n). 

Multiplying k copies of the class of transpositions, we obtain 



( 11 ) 


k times 


where clearly hk }l i{n) = /ifc )7r (n) for any it £ C^n). In other words, (ifc l7r (n) depends on n only up to 
conjugacy class. That this also holds for Weingarten numbers is not so obvious. To see that Weingarten 
numbers are central, we consider the enumeration of strictly primitive factorizations, i.e. factorizations 


7T = (si il) O (s 2 f 2 ) o ■ ■ • O (s k tk) 


( 12 ) 


such that t \ < t '2 < • • • < tk- One may show by a direct combinatorial argument that any permutation 
7 r £ admits a unique strictly primitive factorization, and that this unique factorization has length // 1 . 
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This combinatorial fact can be written algebraically as follows. Consider the Jucys-Murphy elements 
Ji, J 2 , ■ ■ ■, Jt, • ■ ■ £ C[5(oo)] defined by 

Jt = $>*)• (13) 

S<t 

Let denote the alphabet {{ Ji, J 2 , ■ ■., J n , 0,0,... }}. Then 


E ^{length k strictly primitive factorizations of 7r}7r 

7rG5(n) 

h k (^-n) = E ^{length k primitive factorizations of 

7t G5(n) 


(14) 


where 


— y 1 1 ... Xt k 

t 1 <t 2 <--<t k 

hk= E x tl x t 2 ...x tk 


(15) 


are the elementary and complete symmetric functions in commuting variables The fact that 

each 7 r G CL admits a unique strictly primitive factorization, and that this factorization has length 7 r, 
translates into the identity 

efc(S n ) = E G v( n ) e 2 {n), (16) 

jW=fe 

which was first obtained by Jucys HI (see also SI)- On the other hand, the algebra A of symmetric 
functions is precisely the polynomial algebra A = C[ei, e 2 ,... ] in the elementary symmetric functions, 
so we conclude that the substitution / 1 —»• /(H„) defines a specialization A —> Z[n) from the algebra of 
symmetric functions to the class algebra. In particular, /ifc(£ n ) G Z(ri), and we can write 


ftfc(Sn) = E w k,n(n)C li (n), 


(17) 


where u>k,g,(n) = Wk^{n) for any 7r G C M (n). 

2.2 Character theory 

Since any permutation is either even or odd, the Hurwitz and Weingarten numbers ,,(«), Wk,i, (ji) can 
only be non-zero for k of the form k = |// | + 2g for integer g > 0. We thus introduce the notation 

hg,n{n) := ll|^|+2 

WgAn) ■= ^| M [+2 gA n )- 

In particular, Theorem[I]reads 

(~l)MN n+ M W«rr = En^ 

s>0 


( 19 ) 
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where n < N and o o r 1 £ Using the character theory of S(ri), Jackson (7) and Shapiro- 

Shapiro-Vainshtein ED obtained the remarkable formula 

1 n ~ 1 / -| \ / / \ \ n— l+2p 

/i s ,(n-i)(n) = ^5Z(- 1 ) J ( n i (20) 

for the number of factorizations of a full cycle (i.e. an element of C( T j_ 1 )(n)) into n—l+2g transpositions. 
Here we will explain how properties of Jucys-Murphy elements in irreducible representations of C[S(n)} 
may be used to obtain an analogous formula for the Weingarten number w g ^ n _i){n). 

Our point of departure is the remarkable expansion 

/( 2 „) = E ( 21 ) 

obtained by Jucys ©, of the symmetric function / £ A evaluated at S„ in terms of the characters 

X A :=^ X "(^(n))C» (22) 

of the irreducible (complex, finite-dimensional) representations of C[<S(?i)]. Here A\ denotes the alphabet 
of contents of the Young diagram A, and H\ is the product of its hook-lengths. This can be viewed as 
an analogue of the formula of Burnside which expresses the connection coefficients of Z ( n ) in terms of 
irreducible characters. 

Consider the ordinary generating function 

$0; n) = hk(Z n )z k , (23) 

k> o 


which is an element of the algebra Z(n)[[ 2 ]] of formal power series in one indeterminate z over the class 
algebra Z(n). Plugging in the character expansion 


MH„) = Y. 


Ahn 


and changing order of summation, we obtain 

®(z;n) = E 




X 


£(□)*)’ 

where c(D) denotes the content of a cell □ £ A and we have made us of the generating function 


E hk(X!,X 2 , ■ ■ ■ )z k = — 


1 


1 — XsZ 


k> 0 i>l 

for the complete symmetric functions. Thus we obtain the formula 

X A (CV(™)) 


Ahn 


ineA 


(l-c(D)z) 


(24) 


(25) 


(26) 


( 27 ) 
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for the ordinary generating function 


$^( 2 ;; n) = Y w k ^(n)z k 


(28) 


fc >0 


of Weingarten numbers. Note that by Theorem[I] this corresponds to the character expansion 


W - = E 


X X {C ll {n)) 


Ah n 


HxUuex( N + c ( a )) 


(29) 


(where n < N and a o t ~ 1 £ C^(n)), which is well known in the physics literature and was first 
rigorously obtained in El by a different argument. 

Up until this point, the partition // has been generic, but now we restrict to the special case fi = (n— 1), 
the class of a full cycle in S(n). A classical result from representation theory informs us that the trace of 
C , ( Tl _i)(n) in an irreducible representation can only be non-zero in “hook” representations: 


X X {C( n - 1 ) 00 ) 


(-l) r , if A = (n — r, l r ) 

0 , otherwise 


(30) 


Now, the content alphabet of a hook diagram may be obtained immediately, 
A 

so that 


n-r, 1 -) = {0,1, ..., n - r - 1} U {-1, ..., -r}. 

(-1 r 


n— 1 


$(„_!) {z; n) = Y 


ho H (n- r ,in nr=i (! - iz ) rij-h 1 +j-) 

For example, if n = 4, this is a rational function of the form 

const. const. 

^ ’ (1 — z)(l — 2z)(l — 3 z) (1 — z)( 1 — 2z)(l + z ) 

const. const. 

(1 — z)(l + z)( 1 + 2z) (1 + z)(l + 2z)(l + 3s) 

Thus, as an irreducible rational function, $(„_!) (z; n) has the form 

— 1 i 

< f > ( n _i)(z; n) = — ~* =0 Cl4/ - 

nirh 1 -^ 2 ) 

where cq,..., c„_ 1 £ C are some constants to be determined momentarily. 

Before finding the above coefficients, let us consider the generating function 

1 


n’uu-i 2 «) 


(31) 


(32) 


(33) 


(34) 


( 35 ) 
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The coefficients in this generating function are complete symmetric functions evaluated on the alphabet 
{l 2 , ... ,n 2 } of square integers. Reason dictates that they ought to be close relatives of the Stirling 
numbers 

S(n + g,n)=h g {l,...,ri). (36) 

The Stirling number S(a, b) has the following combinatorial interpretation: it counts the number of parti¬ 
tions 

{1, • • ■, a} = Vi U • ■ ■ U V(, (37) 

of an o-element set into b disjoint non-empty subsets. Stirling numbers are given by the explicit formula 


S(o,6) = E(" — 


3= 0 


j)V 


(38) 


The numbers 

T(n + g,n) =h g (l 2 ,...,n 2 ) (39) 

are known as central factorial numbers. The central factorial numbers were studied classically by Carlitz 
and Riordan, see 1lTj Exercise 5.8] for references. They have the following combinatorial interpretation: 
T(a, b) counts the number of partitions 

{l,l , ,...,a,a , } = ViU---UV r 6 (40) 

of a set of a marked and a unmarked points into b disjoint non-empty subsets such thaj^J for each block 
Vj, if i is the least integer such that either i or i' appears in Vj , then { i , *'} C V :j . Central factorial numbers 
are given by the explicit formula 


T(a,b) = 2^(-l) 


b ~j. 


■2a 


3=0 


(b-j)\{b + j)V 


(41) 


Now let us determine the unknown constants cq, ..., c„_i . By the above discussion, the generating 
function •!>(„_!) (z; n) has the form 


$(„_!) (z;n) = (c 0 + c\z -1-c„_i z n x ) ^ T(n - 1 + g, n - l)z 2a . (42) 

g> o 


On the other hand. 


$(„_!) (z;n) = Wk,( n -i){n)z k 
k> o 

= E<("-i )^)z n ~ 1+29 ^ 

g> o 

= w 0 ,{n-i)(,n)z n ~ 1 + Wl,(n-l)(n)z n+1 + .... 


® Balint Virag (personal communication) gave a colourful description of this condition, which is actually quite a useful mnemonic: 
“the most important guy gets to bring his wife.” 
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Consequently, we must have co = • ■ ■ = c n _2 = 0, c n _i = Wo,(n-i)( n )i the number of primitive factor¬ 
izations of the cyclic permutation £[ 1 , n\ = (1 2 ... n) into the minimal number n — 1 of transpositions. 
It is not difficult to show (see l5l fl0l ) bijectively that the number of minimal primitive factorizations of 
the cycle £[l,n] is the Catalan number Cat n _i = -( n “ 2 )- In fact, a stronger result from fTOl asserts 
that the number wq ^(n) of minimal primitive factorizations of an arbitrary permutation of reduced cycle 
type fi is a product of Catalan numbers, 

Am) 

wo,,j.{n) = JJ Cat Mi , (44) 

*=1 


so that the function 


7 r 1 —t #{minimal primitive factorizations of n} (45) 

is a central multiplicative function on S( 00 ) (note that, via Theorem [I] this result corresponds to the 
first-order estimate 

Am) 

=n Cat Mi +0 

i=l 

where eror -1 £ C M (n)). Thus we have proved the following analogue of ( |20| ) for primitive factorizations. 

Theorem 2 For any g > 0, the number of primitive factorizations of a full cycle from S(n) into n — 1 + 2g 
transpositions is 



w g ,(n-i){n) = Cat n _i -T(n -l + g,n-l), 

where Ti^a, b) denotes the Carlitz-Riordan central factorial number. Equivalently, we have the generating 
function 




Cat„i z r ‘ 


(1 — l 2 z 2 ) ... (1 — (n — l) 2 z 2 )' 

Via Theorem[l] Theorem|2]corresponds to the exact integration formula 

(—l )" -1 Cat„! 


Wr7T N(N 2 — l 2 )... (TV 2 — (n — l) 2 ) : 


(JOT ^ £ C,j,(n), 


which was first stated by Collins in J2)- 


(47) 


3 Conclusion 


We have discussed the close relationship between the problem of computing polynomial integrals over 
unitary groups and the enumeration of primitive factorizations of permutations. In particular, the problem 
was completely solved for full cycles, and the central factorial numbers of Carlitz and Riordan made a 
surprising appearance and were given a new combinatorial interpretation. It seems that Hurwitz numbers 
and Weingarten numbers are remarkably similar in character. For example, writing ( |20| > and Theorem [2] 
in terms of exponential generating functions yields 


hg,( n -i)(ri) = n n 2 n 2g 


g,(n—l)(jl) — Cat n _i 


n — 1 + 2g 
n — 1 

2n — 2 + 2g 
2n — 2 


-2 g 


(25)!. 

- z l9 

.Mi 


/ sinhz /2 
/ sinhz/ 2 N 2n ~ 2 

\~^J2 


( 48 ) 












Unitary Matrix Integrals, Primitive Factorizations, and Jucys-Murphy Elements 


411 


On one hand, the multiplicative form of Theorem [2] suggests the existence of an underlying bijective ex¬ 
planation, and on the other computer calculations performed by Valentin Feray (personal communication) 
suggest that such a bijection could be very complex. Further similarities between Hurwitz numbers and 
Weingarten numbers are the subject of work in progress 0. Let us finish by pointing out that the first au¬ 
thor has extended many of the results presented here to the setting of polynomial integrals over orthogonal 
groups J9]- 
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